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Abstract 

We examine perturbations of eigenvalues and resonances for a class of multi-channel quantum 
mechanical model-Hamiltonians describing a particle interacting with a localized spin in dimension 
d = 1,2,3. We consider unperturbed Hamiltonians showing eigenvalues and resonances at the 
threshold of the continuous spectrum and we analyze the effect of various type of perturbations 
on the spectral singularities. We provide algorithms to obtain convergent series expansions for the 
coordinates of the singularities. 
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1 Introduction 

An extensive recent literature presenting different rigorous approaches to the analysis of perturbations 
of energy eigenvalues embedded in the continuous part of the spectrum of Schrodinger operators is 
now available (see e.g. [6], [7], [8], [9], p3], [H] and references therein). 

Less is known about the case of eigenvalues embedded at the threshold of the continuous spectrum. 
Mainly because dilation-analyticity techniques loose their effectiveness when applied to the study of 
perturbations of bound states or resonances at a threshold, results on this particular case are rare. At 
the best of our knowledge the only recent work on this topic was done by A. Jensen and G. Nenciu 

They consider the case of a Schrodinger operator 

H{e) = -A + V + eW = H + eW 

on L 2 (M n ) with n odd. The unperturbed Hamiltonian H is assumed to have a non degenerate eigen- 
value at zero energy and suitable hypotheses on V guarantee that the essential spectrum of H is 



purely absolutely continuous and fills the half line [0,oo). The self-adjoint operator W is assumed to 
have strictly positive expectation value in the eigenvector of H at zero energy, in such a way that the 
singularity is "pushed up" by the perturbation. 

Under some technical assumptions on the properties of {H — z)^ 1 for complex z, close to the origin, 
the authors prove that, as effect of the perturbation, the zero eigenvalue develops into a resonance. 
It is known that the perturbation drives an eigenvalue strictly embedded in the continuum spectrum 
of H to move into a singularity of the resolvent in a non-physical Riemann sheet. The imaginary part 
of the singular point position grows in that case like e 2 for small s. On the other end, in the case 
of an eigenvalue at threshold, Jensen and Nenciu find a behavior like e 2+u ^ 2 , with v integer v ^ — 1. 
This implies that the lifetime of the corresponding resonances shows a universal dependence on the 
perturbation strength in the former case whereas in the threshold case it can be larger for v 1, 
or smaller in the case v = — 1. In several examples of one and two channel Schrodinger operators in 
dimension 1 and 3 the authors detail the computation of the leading order in e of the lifetime of the 
resonance. They show explicit cases where v 1 and others where v = — 1. 

In the following, in order to investigate formation of resonances by perturbation of threshold eigenval- 
ues we will make use of Hamiltonians characterizing different dynamical models for a quantum particle 
moving in an array of localized spins (see [3] for details). By changing their geometrical and dynamical 
parameters the spectral structure of the Hamiltonians can be adapted to show isolated or embedded 
eigenvalues as well as eigenvalues and resonances at any threshold of the continuous spectrum. In 
this paper we analyze only two channel Hamiltonians for one particle and one spin with eigenvalues 
and/or resonances in the upper channel at the continuum threshold. 

Inasmuch as we examine specific models our analysis lacks some generality. However we want to 
point out that zero-range interaction Hamiltonians are particularly versatile models reproducing all 
the qualitative dynamical features typical of short range potentials (see e.g. [2], [5]). In particular one 
and many channel point interaction Hamiltonians share the property of having a resolvent expansion 
around the origin of the type assumed by Jensen and Nenciu to prove their results. Conversely any 
Schrodinger operator having such resolvent expansion around the origin can be approximated by point 
interaction Hamiltonians. Moreover the high degree of computability typical of the models we discuss 
here allows us to consider also two dimensional cases and write down explicitly convergent series 
expansions for the coordinates of resonances and eigenvalues. 

In this paper we will only consider multichannel Hamiltonians in order to have a more direct compari- 
son with the results obtained for embedded eigenvalues in [I]. One channel Hamiltonians, possibly with 
multiple point-scatterers, also show very rich spectral configurations with resonances and eigenvalues 
at some continuous threshold. We plan to examine those cases in further work. 

We will not state here our results in a time-dependent framework, limiting ourselves to the so called 
spectral form of the Fermi Golden Rule. The investigation of the "survival probability" of the resonant 
state is a fundamental step in order to investigate reality and time range of validity of the expected 
exponential behavior on which the very notion of lifetime relies. We mention that the complete 
knowledge of the Hamiltonian generalized eigenfunctions allows in our case a very detailed analysis of 
the time evolution of the resonant state as it was done in the case of embedded eigenvalues ([!]). 
The paper is organized as follows. In section[2]we introduce notation and basic definitions. In section[3] 
we state and prove our results. Within this section we split in subsections the analysis of different kinds 
of perturbations. A final section consists of a summary of results together with further comments. 

2 Basic definitions and results 

For d = 1,2,3, we consider the Hilbert space := L 2 (R d ) L 2 (R d ). We denote by f the generic 
(column) vector in Jif: 





2 



is the state space of a quantum particle in M rf in presence of a two-level quantum system (a spin) 
localized at the origin, i/jq ( resp. ipi) represents the state of the particle in the channel where the 
spin (better a particular component of the vector spin operator) has value — 1 (resp. +1). In what 
follows we consider Hamiltonians in ^ belonging to the family of the self-adjoint extensions of the 
symmetric operator S defined by D(S) := Cg°(K d \{0}) C§°(K d \{0}), Stf := (-Aipo, (-A + l)V>i)- 
We do not detail here how to characterize the whole family of self-adjoint extensions of S, as this was 
done, in a more general setting and with slightly different notation, elsewhere (see [3] for d = 1, 3, and 
[3] for d = 2). According to the definition of S the 0-channel (resp. the 1-channel) will be referred to 
as the lower (resp. the upper) channel. 

For z G C\M + , we denote with G z (x) the fundamental solution in M. d of Helmoltz's equation: (—A — 
z)G z = 5; explicitly 



G z {x) 



e ly/z\x\ 

i 1=- 

e iyfz\x\ 

Att\x\ 



d = l 

d = 2 
d = 3 



with z G 



; Im (v£) > 



(2.1) 



Here Hq 1 ^ (77) is the zero-th Bessel function of third kind (also known as Hankel function of the first 

kind), see, e.g., [I]. We recall that Hq 1 ^ (77) tends to zero as \rj\ — > 00 for Im 7/ > and that it has a 
logarithmic singularity in zero 



H, 



(i) 



'jl) 



2*1^5 + 1 + ^ + 0(ln(7y)7y 2 ) 

7T Z 7T 



where 7 is the Euler's constant (7 ~ 0.577). 

Notice that, for Im (\fz) > 0, G z (x) G L 2 (M. d ) for d = 1, 2, 3, a property which does not hold in higher 
dimensions. This is a crucial feature in the definition of our model-Hamiltonians and it is the reason 
why point interaction Hamiltonians are trivial in dimensions bigger than three. 
We denote by Hq the Hamiltonian in Jt? given in the following 

Definition 1. Let 9q G E and 
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d 
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ln(2) - 7 



~ 54 . 



(2.2) 



Hq : D{Hq) C J^ 7 — )• ffl 1 is the self-adjoint operator: 



D(H ) :-- 



* = 



00 

i 4 ; i 



G ^ 



0o = <P Z + (?oG 2 , 0i = 0? + giG 2 - 1 ; 0§ , 0f G -ff^(IR d ) ; * G 



r2 /TTD<i\ 



qo — %fo , 91 — 0i /1 ; 

/i=^(0) 



fi = lim 

3 \x\-+0 

fi = lim 

3 |*|-+o 



^■(x) + |^ln(|x|) 



V>j(x) 



2tt 
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j = 0, 1 , d = 2 ; 

J = 0, 1 , d = 3 
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For any choice of 6q E M, Hq belongs to a sub-family of the self-adjoint extensions of the operator S 
defined above. Each Hamiltonian of the sub-family generates a dynamics where the spin evolution is 
not affected by the particle ( |4] ). The two channel Hamiltonian Hq is often formally written as 



-A + a 8 











-A + l + ai6, 



to stress that it acts as S on D(S). To match our notation with the one used in the monograph [2] 
on point interactions one should take into account the following correspondence rules: ay = — 9j for 
d = 1; ay = 1/0 j for d = 2, 3; j = 0, 1. 

The spectrum of Hq can be obtained directly from the spectrum of the operator "—A + (XjS" ( see 
[2]). The main results are collected in the following 

Proposition 1. For d = 1,2,3, the essential spectrum of Hq fills the positive real line, a ess {Ho) = 
[0,oo) and E o~ p {Hq), the point spectrum of Hq. The bound state (of unit norm) corresponding to 
the zero energy eigenvalue is $ = (0, <\)\ ) with 



i\x\ 



V2 



-kl 



47r|x| 



d 
d 

d 



Moreover: 

- Ford=l and 6 > 0, a p (H ) = {-6^/4,0}. While for 0q ^ 0, a p (H ) = {0}. 

- Ford = 2 and 6 G R, 6 ^ 0, a p (H ) = {-e 4 ^ 1 /"-!/^] , 0} (where the constant a was defined in 



equation (2.2)). While for 6q 
- For d = 3 and 9q < 0, a p (Ho) 



0, a p (H ) = {0}. 

= {-(4^ ) 2 ,0}. While for 6 > 0, a p (H ) = {0}. 



As noticed before Hq describes a two independent channel system. In each channel the particle "feels" 
a point interaction placed in the origin whose strength may depend on the channel (equivalently on 
the spin state). Among all the self-adjoint extensions of S one can find a large class of Hamiltonians 
coupling the two channels. To the aim of examining the behavior of the eigenvalue of Hq at the 
threshold of the essential spectrum, when the lower and upper channels are weakly coupled, we choose 
a suitable Hamiltonian, H e , belonging to that class and close, in a sense that will be made precise in 
the following, to Hq. 



Definition 2. Let us take 6q, b, c G E and let e > 0. For d 

( 2 + ce d=l 



a + cs 
— An + ce 



d 
d 



1, 2, 3 we set 



2vr 



ln(2) - 7 



H £ : D{H £ ) C — > is the self-adjoint operator: 



D(H e ) :-- 



* = 



4>o = <P z q + qoG z , Vi = <f>{ + qiG z 



z E 



90 : 

/i 
/i 

/i 



^o/o + te/i , qi = befo + e\h ■ 



lim 

|a|->0 

lim 

|aj|->0 



^(x) + ^ln(|x|) 



0,1, d 
0,1, d 



47r|x| 



i = 0, 1 , d = 3 
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-Acf> z + zqoG z 
+ + z qi G z - 1 ' ' 



G D{H e ) 



For all z G C\R we denote by R £ (z) := (il £ — z) -1 the resolvent of H £ . An explicit formula for R e (z) 
can be obtained by using the theory of self-adjoint extensions of symmetric operators (see f3] and [1]). 
We summarize the result in the following formula 



Re(z) 





(-A + l-z)- 1 
T £A1 (z)(G z ,-)G z 

K T £ , 21 {z)(G z ,-)G z - 1 
where the function G z was defined in equation (|2.1l) and 



De{z) 




T £il2 (z)(G z ~ l ,-)G z 
T £ , 22 {z){G z -\-)G z - 1 



(2.3) 



' b 2 e 2 - (6» + 2iyfz) [2(1 + Wz - 1) + ce] d = 1 

{a + ce+ [6 (a + ce) -b 2 e 2 ]g{z)}{9 + [6 (a + ce) - b 2 s 2 ]g{z - I)} - b 2 e 2 d = 2 



i + t 



2-1 



where (7(2) := [ln(y^) — z7r/2] /(27r) — 1/a and the matrix elements Y £ ^j(z) read: 
for d=l 

T £iU (z) = -2iy/z [-b 2 e 2 + 6> (2 + 2iVz - 1 + ce)] 



d = 3 
(2.4) 



for d=2 



r £) i2 (2) = r £)2 i (z) = 4be V^T 

r £ ,22(^) = -2iVz~^l [-b 2 e 2 + (2 + ce)(2»v^ + O )] , 

r e ,n(2) = {^o + [Oo (a + ce) - b 2 e 2 ]g{z - l)}[9 (a + ce) - b 2 e 2 } 

FeM z ) = r e,2i(^) = be[9 (a + ce) - b 2 e 2 } 

F £ , 22 {z) = {a + ce+ [9 (a + ce) - b 2 e 2 }g(z)}{9o{a + ce) -b 2 e 2 } 



(2.5) 



(2.6) 



for d=3 



T £ , n (z) = 9 
^e,n(z) = T £t21 (z) = be 
F £ M Z ) = (-4vr + ce) 



CE 



6 z e 



47TZ 



(2.7) 



1 • ^° A" 



47U 



For d = 1,2,3 the explicit form of the resolvent of -fTo, Ro(z) '■= (Hq — z) , can be obtained from 
formulas (2.3) - (2.7) by setting e = 0. 

Remark 1. For 60 = and d = 1,2 the resolvent #0(2) has a singularity in z = on both channels. 
In the upper channel there is a polar singularity corresponding to the eigenvalue. In the lower channel 
there is a singularity of order z~ 1 ^ 2 for d = 1 and lnz for d = 2 respectively. This is easily checked 
analyzing the behavior around z = of (—A — z)^ 1 (see, e.g., [12J). A precise statement, obtained 
examining the integral kernel G z (x — y) of (—A — z)~ x , gives the following expansion 



(-A-z)- 1 



2jz 
1 

47T 
10(1) 



+ 0(1) 



In* + 0(1) 



d= 1 

d = 2 
d = 3 
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where 0(1) denotes an operator on some suitable weighted L 2 space, whose norm remains bounded 
uniformly in z. A possible choice for the weighted space is for example L 2 (R d , (1 + \x\)~ s dx) for some 
s large enough. 

Finally we notice that H £ is a small perturbation of Hq in the resolvent sense, i.e., \/z £ C\R there 
exists £q such that for all < e < £o 

\\R £ (z) - Ro(z)\\b(j?,j?) ^eC 

where C is a positive constant independent on e and || • \fs{& &\ is the operator norm in the vector 
space B(J4f, Jff) of bounded linear operators on 



3 Results 

In this section we analyze the spectral structure of H £ to examine the effect of the coupling between 
the lower and the upper channels on the spectrum of the Hamiltonian Hq and in particular on its zero 
energy eigenvalue. We denote by a p (H £ ), a ess (H £ ) and a ac {H £ ) the point, essential and absolutely 
continuous spectrum of H £ respectively. 

It will be clear from the proofs that value and sign of the parameter b in definition [2] do not affect our 
results in any substantial way. For this reason we set 6=1. 

We use sometimes the phrase "for e small enough ... " as a short version of "there exists £q such that 
for all < e < e ... ". 



3.1 Positive perturbations, c < 0. 

In this section we study the behavior of the threshold eigenvalue when the parameter c in the Hamil- 
tonian H £ is negative. This choice corresponds to a positive perturbation of the Hamiltonian Hq in 
the sense that for small e the threshold eigenvalue is pushed inside the continuum as a consequence 
of the perturbative term ce in B\. 

In order to make this statement more precise let us consider the case b = in definition [2] The 
Hamiltonian H^ =0 is a perturbation (in resolvent sense) of the Hamiltonian Hq for which the channels 
and 1 are not coupled. For all e small enough the Hamiltonian H^ =0 has an eigenvalue E^ =0 in a 
ball of radius e around the origin, and 

ae2 2 2 

E b=o = 1 _^ = _ c£ _^_ d=1 (31) 

fc 4 4 

E b=o = j _ e 4*(i/«-i/flf) = + , 2) d = 2 (32) 

a 1 



ET u = l-^ = -^ + 0{e') d = 3. (3.3) 



Looking at H £ as perturbation of H^ =0 , we expect that the zero energy eigenvalue of Hq will be driven 
in a resonance as it happens for embedded eigenvalues. 
We analyze first the cases d = 1, 2 and 9q ^ 0. 

Theorem 1. Let d = 1, 2 and assume that c < and 0q ^ 0. Then there exists Eq > such that for 
all < £ < Sq : 

the essential spectrum of H £ fills the positive real line and is only absolutely continuous, 

a ess (H £ ) = a ac {H £ ) = [0, +oo) ; (3.4) 

there exists a positive constant C such that the Hamiltonian H £ has no isolated eigenvalues in (— C, 0); 
the analytic continuation of the resolvent R £ (z) through the real axis from the semi-plane Im z > 



6 



has a simple pole (resonance) in z = E £ where Im (E]l) < and 
1_ c 2 \ o .2, 



El 



c e + 



er 

7 



el + 0(e 4 ) 



Re (££) 



4-7T I c| 



e + 0(e 2 



Im {El) 



16vr d 



\a9o\ 2 | lne| 



+ o((e/\lne\f 



d=l 



d = 2. 



(3.5) 
(3.6) 



Remark 2. Theorem 1 does not characterize the entire spectral structure of the Hamiltonian H £ . 
Statements only concern spectral singularities in a small region of the complex plane around the origin. 
In particular for 9q > in d = 1 and for 0q 7^ in d = 2 there are negative eigenvalues close the 
corresponding eigenvalues of H^ =0 . The relative singularities of the resolvent of H £ are bounded away 
from the origin unless #0 is very small. 

Looking at the dependence on #0 — > of the coordinates of the resonances in d = 1,2 given in (3.5) 
and (3.6) one realizes that the case 6q = has to be treated independently. 

Proof. We first consider the case d = 1. We notice that the singularity in z = of order z~ 1 / 2 , in the 



term (—A — z) 1 in formula (2.3) (see the remark [I]), is canceled by an opposite singularity arising 
from the coefficient T £ ^i/D £ . In fact an explicit calculation gives 



r e ,ii(.s)(GV)G' 



2^~z 



+ 0(1) 



for all e > 0, where the equality has be intended in some weighted L 2 (IR) space, where the weight 
depends on the number of terms of the expansion one considers (see [12] for details). A similar remark 
holds true for the singularity in z = 1 of order (z — l)" 1 / 2 in the upper channel. This singularity, 
arising from the term (—A + 1 — z)^ 1 in the formula (2.3), is compensated by an opposite singularity 
in the term T £ ^2/D £ . Then the singularities of the resolvent on the real axes coincide with the zeros 
of the function D e {z) in equation (2.3). 

We prove first statement (3.4), showing that H £ has no embedded eigenvalues or eigenvalues at the 
threshold. Let us set z = A > 0. Trivially D e (l) 7^ and D £ (0) ^ 0. For A G [0, 1) a direct calculation 
shows that equations Im [Z? e (A)] = and Re [Z) e (A)] = are not compatible. We deduce that there 
are no solutions to the equation D e (A) = for A G [0, 1]. For A > 1, taking real and imaginary part 
of the equation D(X) = we get 

(2 + ce)V\ = -OqVX - 1 



and 



£ 2 + 4\/A VA - 1 = #o(2 + ce) 



(3.7) 



(3.8) 



For 6*o > equation (3.7) has no solutions in (1, +00) and for 6>o < equation (3.8) has no solutions 
in (1, +00). 

Next we prove that there are no isolated eigenvalues in some suitable neighborhood of z = 0. Let us 
set z = —A, then equation D £ {—\) = gives 

-6 + 2\/A = . (3.9) 

2(1 - x/A+T) - \c\£ K ' 

For A G (0, +00) the right side of the equation is a strictly positive and decreasing function which 
equals e/\c\ for A = 0. While the left hand side of the equation is a strictly increasing function which 
equals —6q for A = 0. Then for e small enough there are no solutions of (3.9) when #0 < 0. For 9q > 
the l.h.s. has a zero in A = Oq/A and it is negative in (0, # 2 /4) and positive in (0q/4, 00). Then for 
#0 > there is only one solution to equation (3.9), say Ao, e , and Ao )£ > #o/4. It follows that for #0 / 
there are no isolated eigenvalues in (— 0g/4, 0). 

To find a solution to the equation D £ (z) = we make use of the following recursive procedure. We 
first notice that the equation D £ (z) = can be written as 



iy/z 



1 



-1 



ce 
Y 



2(fl + 2iy^) ' 
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We look then for a fixed point of the recurrence relation 
z<® = 



z 



(fc+i) 



1 



1 



c e 



e 



1 2 



2 2(6> + 2i\/^W) 



fc = 0,1,2,. 



(3.10) 
(3.11) 



To prove convergence of the sequence z^ we proceed by induction. Assume that Ce, then 

equation (3.11) implies that for e small enough \z^\ < Ce. Since z^ = |c|e + 0{e 2 ) then \z^'\ ^ Ce 
for all k. Moreover let us set z^ = \c\e(l + w^), k = 1, 2, 3, .... Trivially |w;W| ^ Ce and by equation 
(Hit 



1 

|c|e 



2 1 



ce 



2 / V2((9o + 2i v / Hev / l + w (fc) ) 2(0 O + 2i v / j^jev / l + w^- 1 )) 



1 

Icle 



(3.12) 



.4(6*0 + 2iy / |c|e\/l + ^( fc )) 2 4(^o + 2i v ^e v / l + w^" 1 )) 
Since the function (#o + 2i^\c\e\/l + w)^ 1 is analytic for w in a ball of radius e around the origin 

1 1 <^ Ce^lw^ - w^\ . 



2((9 + 2i v / |^|ev / l + w( fc )) 2(9 + 2i^/\^e\J\ + w^-V) 



The second term in the r.h.s. of equation ( |3.12 ) can be treated in a similar way Then \w^ k+1 ^ — ^ 
Ce z l 2 \w^ -w^' 1 ^ for all k = 2,3, ... which in turns implies that - z^\ < Ce 3 / 2 |,z( fc ) -z^ 1 )] 

for all k = 2, 3, ...; the sequence {z^} converges in a ball of radius e and 



,(2) 



c e + 



e2 +0(e 4 ). 



Let us now consider the case d = 2. We notice that, similarly to what happens in the case d = 1, the 
logarithmic singularities in z = and z = 1 due to (—A — z)" 1 and (—A + 1 — z)~ l in equation (2.3) 
(see the remark [I]) are compensated by opposite singularities arising from the coefficients T £ ^u/D e 
and T Et 22/D e . Then the singularities of the resolvent on the real axes coincide with the zeros of the 
function D e (z). 

First we prove the statement (3.4). Let us set z = A > and analyze the equation D £ (\) = which 
can be written as 



o + ce + (6> (a + ce) - e 2 ) ( 



/ In A 



1 

a 



e + (9 (a + ce)-e 2 )(^^ 



4-7T a 47 

If A > 1, taking the imaginary part of both sides of the last equality one gets after few manipulations 

-21 2 



-{00 + [9 (a + ce) - e 2 ] ( ln(A - l)/(4vr) - 1/a) }' 



[6> (a + ce) - e 2 
16" 



+ e z 



A > 1 



which obviously has no solutions. For < A < 1 one can see that equations Im(Z) £ (A)) = and 
Re(-D e (A)) = are not compatible. Moreover, as we already noticed, the resolvent R £ (z) does not 
have singularities in z = and z = 1. Then the equation D e (\) = has no solutions in (0, oo). 
As a second step we look for isolated eigenvalues of H e in some suitable interval (— C, 0). For A > 
D £ {— A) is real and equation D e {— A) = can be rearranged as 



{6 + [e (a + ce)-e 2 }g(-X-l)} 



a + ce + [9o(a + ce) — e 2 ]g(— A) 



A > 0. 



(3.13) 
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Let us consider first the case 9q > 0. The l.h.s. of equation ( |3. 13 ) is continuous and strictly increasing, 
moreover it equals 6o\c\e/a + e 2 /a for A = 0. The function on the r.h.s. of the equation has a vertical 
asymptote in A flj£ = exp [47r((0o/a— l)(a — |c|£)— £ 2 /a) / (0o(a — \c\e) — £ 2 )] . It is continuous and strictly 
decreasing in (0, A ae ) and (A a)£ ,oo). Finally, in (A aj£ ,oo), it is strictly positive and its limits in A = 
and A = oo equal zero. The above discussion makes clear that equation (3.13) has always one solution 

exp[47r(l/a-l/0 o )] 
A similar 



in (X a ,£, oo) and that such solution is contained in a ball of radius e around A aj o 
which is the eigenvalue of the Hamiltonian Hq. For c < there are no solutions in (0, A a)£ ) 
result can be proved by an analogous argument if one assumes 9q < 0. We have then shown that for 
c < and #0 7^ 0, H £ has no isolated eigenvalues in (— C, 0) with < C < A a)£ . 

Let us now prove that for c < and 0o 7^ the equation D e (z) = has a solution with negative 
imaginary part in a neighborhood of z = 0. Again we make use of a recursive procedure. The 
equation D £ (z) = can be rearranged as follows 



9{* ~ 1) 



00 



a — c £ 



)-- 2 ]{ 



a — c £ 



+ [0o ( 



a — c £ 



a — c £ 



(3.14) 



Then for all k = 0, 1, 2, . . . we define the sequence {*(*)} by 

z<® = 



ln(l - z {k+ ^) = r(zW) k = 0,1,2 



(3.15) 



where r{z) is the function 

47r [0 o |c|e + £ 2 ] 



r(z) :- 



+ 



[0o(a — \c\e) — £ 2 ] [0o(a — \c\e) — £ 2 ] {a — \c\e + [0o(a — \c\e) — e 2 ]g(z)} 



Taking the imaginary part of the equation (3.15) one has 



arg(l 



\c\s+ [e (a- \c\e) -e 2 ]g(z( k ))\ 



IT- 



arg(z (/c) )) 



(3.16) 



We remark that the Riemann surface associated to the function arg(z) is made up of an infinite number 
of sheets, each one labeled by an integer number n and characterized by arg(z) G [2mr,2(n + 
We call "first" (or "physical") Riemann sheet the one corresponding to n = 0, i.e., arg(z) G [0, 27r). 
If arg(z) G [— 2tt, 0) we say that z is in the "unphysical" Riemann sheet corresponding to n = — 1. 
In each sheet of the Riemann surface arg(-^/z) = arg(z)/2. If z is in the "unphysical" Riemann sheet 
corresponding to n = —1, then Im (y/z) ^ 0. 

In equation (3.16) the function avg(z^) has to be intended as the analytic continuation of the argument 
function through the cut [0, 00). Notice that we assumed arg(z^) = 0. 

We proceed by induction. First we prove that there exists £0 such that for all < e < £0, if 
-7T < arg(z( fc )) < and C x e < \z^>\ < C 2 £ then -vr < arg(z( fc+1 )) < and C 3 £ < C 4 £, 
where the constants C3 and C4 do not depend on C\, C2 and £0 = £o(Ci, C2, C3, C4). We first notice 
that if — 7T ^ arg(z( fc )) ^ and C±e ^ \z^\ ^ C 2 £, then there exists £q such that for all < e < Eq 



Aire 2 



C £ 



+ [0o (a 



C £ 



£ 2 Wz( fc ))| 



In el 



(3.17) 



z 



(k+iy 



From equation (3.16) and estimate (3.17) one immediately gets that for e small enough ^ arg(l 



) ^ C£ 2 /|ln£| 2 , consequently —it ^ arg(z( fc+1 )) ^ 0. From the definition of the function r(z) 
and from the inequality (3.17) it follows that |r(,z( fc+1 ))| ^ Ce. Then from the analyticity of the 



exponential function one has 



,(k+l) 



-r(zW) + 0{\r{z^)\ 2 ) 



(3.18) 
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which in turns implies 



(3.19) 



and one can take, e.g., C3 = 2-K\c\e/a 2 and C4 = 87r\c\e/a 2 . 

We prove now that for all k = 1,2, - ^ Ce| ln£|~ 2 |2i( fc ) - z( k ~^\. For fc = 2,3, ... let us 
set = (47r|c|e/a 2 )(l + w^). By equation (|3.18[) it follows that, for any fe = 1, 2,3, 



_ z (*)| 



r(«W)-r(^)| 



where we used the analyticity of the exponential function. Equation (3.19) implies that \w^\ ^ Ce 
for all k = 1,2, .... The function r((4-7r|c|e/a 2 )(l + to)) is analytic for w in a ball of radius e around 
the origin, and 

\r(z (k) ) - r(z( fc - 1 ))| = |r((4vr|c|e/a 2 )(l + w {k) )) - r((4vr|c|e/a 2 )(l + w (fc+1) ))| 

_2 



h,C*)- «,(*-!> 1 



lne| 2 I In e| 2 

Since arg^ 1 )) = and z^ = 47r\c\e/a 2 + C(e 2 ) we have proved that the sequence is convergent 
in a ball of radius e. The lowest order of E T £ given in equation (3.6) can be computed by using formula 
p7l5| ). □ 

In the next theorem we analyze the behavior of the zero energy eigenvalue when #0 = 0, c < for 
d = 1,2. In addition to the eigenvalues in (3.1) and (3.2), H £ =0 has in this case a resonance at zero 
energy in the lower channel. The coupling between channels will turn the resonance into a negative 
eigenvalue. 

Theorem 2. Let d = 1,2, assume that c < and set 6q = 0. Then there exists £0 > such that for 
all < £ < So : 

the essential spectrum of H £ fills the positive real line and is only absolutely continuous, 



a ess {H £ ) = cr ac (H e ) = [0, +00) ; 



(3.20) 



the analytic continuation of the resolvent R £ {z) through the real axis from the semi-plane Im z > 
has a simple pole (resonance) in z = E £ where Im (E £ ) < and 



El 



C £ 



Re (El) 



4tt\c 



£2 + <D(e< 



£ + 0(£ 2 |hi£|) ; lm(E r £ 



7T 

a 



2 £ 2 +o(£ 2 ) 



f/ie Hamiltonian H e has an isolated eigenvalue E £ < and 

,-2 



4|c| 2 
exp 



< E £ < 



47r|c| 4-7T 

— u + — 

£ a 



< < 



d = 1 

d = 2: 



d = 1 
d = 2. 



(3.21) 
(3.22) 

(3.23) 
(3.24) 



Remark 3. For d = 2 and #0 = the Hamiltonian H E has also one eigenvalue which moves to minus 
infinity as e goes to zero. One can interpret the latter eigenvalue as the effect of the perturbation on 
the "eigenvalue" of Hq at —00 in #0 = 0~ (see Proposition 1). The emergence of the two eigenvalues 
(3.23), (3.24), respectively in d = 1 and d = 2, reveals the effect of the channel coupling on the 
resonances at zero energy in the lower channels, whereas the two resonances (3.21), (3.22) appear as 
the effect of perturbing the eigenvalues in the upper channels. 
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Proof. We notice first that also in this case, similarly to the case Oq (see theore m [T|) , the singular- 
ities in z = and z = 1, arising from (—A — .z) 1 anc ^ ( — A + 1 — -z) -1 in equation (2.3), are canceled 
by opposite singularities in the two terms T S) n/D s and r £i 22/-D e respectively. This is true both in 
d = 1 and for d = 2. 

Statement (3.20) can be proved as it was done in theorem [T] A straightforward analysis shows in fact 
that equation D £ (X) = has no solutions for A G [0, oo). 

We discuss now the existence of isolated eigenvalues in d = 1. Let us set z = —A, A > 0. The equation 
D £ (—X) = can be written as 

(3.25) 



2VA 



2(VA + 1-1) + \c\e 

The left hand side in the last equation is a strictly positive, increasing function which equals zero for 
A = 0. The right hand side of the equation is a strictly decreasing, positive function which equals 
s/\c\ for A = 0. It is obvious that there is only one solution A e to the equation (3.25) and that 
< A e < e 2 /(4|c| 2 ) which in turn implies estimate (3.23). 

We prove now the existence of an isolated eigenvalue of H e for d = 2. For A > the equation 
D £ (—X) = can be written as 

m(l + A) _ 1 1 

47r a a + ce - e 2 (ln A/(4tt) - 1/a) 

Let us denote by f r (X) the r.h.s. of the last equation, 



A > 0. 



(3.26) 



ce + e 2 /a-e 2 lnA/(47r) 
a[a + ce + e 2 /a — e 2 lnA/(4-7r) 



The l.h.s. of equation (3.26) is a strictly positive and increasing function which equals zero for A = 0. 
The function f r has a vertical asymptote in A aj£ = exp [Air[a/e 2 — \c\/e+l/a)~\ and is strictly decreasing 
in (0, X a , e ) and (A ai£ ,oo); we notice that X a , e goes to infinity as e goes to zero. Moreover the function 
f r equals 1/a for A = and zero for A = A e = exp [47r( — \c\/e + l/o)l • It is positive for A £ (0, A £ ) 
and negative for A E (A e , X a<e ). 



Then there are two solutions to equation (3.26). One is in the interval (A aj£) oo). The isolated eigen- 
value corresponding to this solution moves towards minus infinity as e goes to zero. The other solution 



is in (0,A £ ). The corresponding eigenvalue satisfies estimate (3.24). 



We investigate now the presence of poles of the resolvent in the "unphysical" Riemann sheet. Let us 
consider first the case d = 1, similarly to what was done in the previous theorem we use a recursive 
procedure. We rewrite the equation D £ (z) = as 



iy z 



1 



Which implies 



1 



1 



|c|e e^ 

+ ir + 4v^ 



|c|e e 
~Y ~ Aiyfz 



(3.27) 



We set z = \c\e(l + w) and define the recursive procedure 
>(°) = 



10 v 



w 



(fc+1) 



,1/2 



2|c|V|c|(l + w( k )) 



\c\e 



.-3/2 



+ 



4iV|c|(l + w;( fc )) 16|c| 2 (l + u>( fc )) 



0,1,2, 



then z^ = \c\e(l + w^) for all k = 0, 1, 2, .... By induction it is easy to prove that for all k = 0, 1, 2, ... 
and for e small enough, \w^\ ^ Ce 1 / 2 which in turns implies \z^\ ^ Ce. Moreover for all k, 
\ w (k+i) _ w {k)\ ^ C£ i/2\ w (k) _ w (k-i)\ xhen \ z (k+l) _ z (fc)| ^ Ce 1 ^ 2 \z^ - zl k -% consequently the 
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sequence {z^} is convergent in a ball of radius Ce and the solution of the equation (3.27) can be 
written as E r e = . By a straightforward computation it is easy to see that 

El = \ c \e ^=ei + ^e 2 + 0(e 3 ) . 

2y^ 4 

In dimension d = 2 the existence of a resonance can be proved in a similar way. The equation D e (z) = 
can be rearranged as 

1 

g( Z - 1) = — — - 

a — \c\e — e^g(z) 

ln(l -z)= An I c l g + £2 g( z ) = . r(z) . (3.28) 
v ; a[\c\e + e 2 g(z) - a] w v ; 

We prove that, if c < 0, there exists a solution (with negative imaginary part) of the last equation in 
a neighborhood of radius of order e 2 \ ln(e)| of the point 4tt\c\e / a 2 . Let us pose z = 47r|c|e(l + w)/a 2 
and let us define the sequence {w^} by 



which implies 



w 



(°) = o 



w (k+1 ^ = -1 + -^-(l - e^ W >) fc = 0, 1, 2, ... 
47r|c|e v ' 

where 

r(to):=r^ — y-(l + w)J. (3.29) 
Taking into account the analyticity of the exponential function one obtains 

w (k+i) = _i_^^f( u; W) + 0(|f( u ,W)|2/ e ) k = 0, 1,2, ... . 
4-7r|c|e 

Noticing that for \w^\ ^ Ce\ hxe\ the following expansion holds true 

r( w W) = - A ^ + 0{e 2 \\ne\). 
By induction one can prove that |u;( fe )| ^ Ce\ lne| for all k = 0, 1, 2, ... . 

Moreover by using again the analyticity of the exponential function and the fact that r(w) is analytic 
in a ball of radius e| lne| one can see that 

\ w (k+i) _ w w\ ^ ^| f(w W) _ r(w^)\ ^ Ce\wW - . 

Then the sequence {w^} converges in a ball of radius e\ lns\ and {z^} converges in a ball of radius 
e 2 | lne| around the point Ai:\c\e/a 2 . Since lm(r(0)) = ire 2 /a 2 + o(e 2 ) (see equations (3.28) and (3.29)) 
to the lowest order we get 

2 2 

Im^ 1 )) = --A— S in(lm(r(0))) + 0(|f(O)| 2 /e) = vre 2 /a 2 + o(e) , 

47r|c|e 47r|c|e 



which implies the second estimate estimate in (3.22). □ 

In the next theorem we analyze the behavior of the zero energy eigenvalue when c < for d = 3. 

Theorem 3. Let d = 3 and assume that c < then there exists Eq > such that for all < e < 
the essential spectrum of H £ fills the positive real line and is only absolutely continuous, 

a ess (H £ ) = a ac (H e ) = [0, +oo) ; (3.30) 
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there exists a positive constant C such that the Hamiltonian H £ has no isolated eigenvalues in (— C, 0); 
the analytic continuation of the resolvent R e {z) through the real axis from the semi-plane Im z > 
has a simple pole (resonance) in z = E £ where Im < and 



El 



\C\£ 



3|c|V 



1 




+ 0(£ C 



(3.31) 



16vr 2 8vr 2 V 2vr' 

Proof. Using the classical results we summarized in remark [T] and the explicit expressions for the 
components of the matrix T e (see equation (2.7)) one can conclude that for d = 3 the singularities of 
the resolvent R £ {z) coincide with the zeros of the function D £ (z) defined in equation (2.4). 
By direct analysis one can see that the equation D £ (X) = has no solutions for A > from which the 
statement (3.30) directly follows. 

The eigenvalues of the Hamiltonian H e are given by the solutions of the equation D £ (—X) 
A > 0. Such equation can be rearranged as 

_2 



for 



1 + ^v/A 

4-7T 



47T 



1 



ce \ 
47rJ 



+ 



(4*0 



rVA 



Vl + A A > 0. 



(3.32) 



For c < and 9q equation (3.32) has no solutions while for c < and 9o < there is only one 
solution in a neighborhood of radius e of the point (4tt/9q) 2 . Then for c < and e small enough there 
are no isolated eigenvalues in some suitable neighborhood of the origin. 

It remains to analyze the existence of solutions of the equation D £ (z) = in the unphysical Riemann 
sheet in a neighborhood of the origin. The equation D £ (z) = can be rearranged in the following way 



ce • 

47T/ 



From which we define the sequence {z^}: z^ 



(4vr) 2 
= 0, 



47T 



y/z — ly/z . 



(3.33) 



ce/(27r) + c 2 e7(47r) 2 2e z 



y/ z (k) _ i v^y 



[1 - ce/(4vr)] 2 



(47T) 2 (] 

( z (fc) _ y z (k) 



1 2 



^ VzW/(47r)) 



(4tt) 



1 



Ct 



/(47r))^(l-»0 o vW/(47r))' 



fc = 0,1,2,..., 



where the formula for z^ k+1 ^ was obtained by solving for z the l.h.s. of equation (3.33). By induction 
one can prove that \z^\ ^ Ce for all k = 0,1,2,.... Moreover by direct computation one can see 
that \z( k+1 ^ — z^\ < Ce 3 / 2 \z( k ^ — z^ k ~^\. Then the series {z^} converges in a ball of radius e and 
E r £ = z^- 00 '; by direct computation one can prove expansion (|3.31 ). □ 



3.2 Pure off-diagonal perturbations, c = 0. 

The case c = marks the boundary between two different behaviors of the threshold eigenvalue under 
perturbation: the evolution towards a proper eigenvalue (c > 0) and the evolution towards a resonance 
(c < 0). The model shows, in this case, peculiar features, strongly depending on the spatial dimension, 
which we will analyze for d = 1, 2, 3 separately. 

If in addition to c = we set 6q = the Hamiltonian Hq shows, for d = 1,2, both a zero energy 

eigenvalue in the upper channel and a zero energy resonance in the lower one. Once more the case 

6o = requires a distinct analysis. 

As it was done in the previous sections we set 6 = 1. 

We study first the case d = 1. 

Theorem 4. Let d = 1 and assume that c = 0. Then there exists eo > such that for all < e < £o ; 
the essential spectrum of H £ fills the positive real line and is only absolutely continuous, 

a ess (H £ ) = a ac {H £ ) = [0, +oo) . (3.34) 

Moreover for e small enough: 
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if@o > the analytic continuation of the resolvent R e (z) through the real axis from the semi-plane 
Im z > has a simple pole (resonance) in z = E r £ and 



q5/2 
? 



e 3 + 0(e 4 



(3.35) 



there exists a positive constant C such that the Hamiltonian H £ has no isolated eigenvalues in 
(-C,0); 



if @0 = the Hamiltonian H £ has an isolated eigenvalue in z = E £ and 



E e 



£ 4 / 3 
"22/3 



+ 0(e 8/3 ) 



(3.36) 



if 9q < i/ie Hamiltonian H £ has an isolated eigenvalue in z = E £ and 

^2 



E K 



+ 0(e a ) . 



(3.37) 



Remark 4. When 6q = the Hamiltonian H £ shows also two singularities at the same distance from 



the origin of the eigenvalue (3.36) and with arg(z = — 5tt/3) and arg(z = —tt/3). 



Proof. Similarly to the case c < 0, see theorems [T] and [2j for any 9q £ M the singularities in z = 
and z = 1, arising from (—A — z) -1 and (—A + 1 — z) , are compensated by the terms T e \\/D e 
and T £) 22/D £ respectively. Then the statement ( 3.34| ) is a consequence of the fact that the equation 
D £ (\) = has no solutions for A > and the singularities of the resolvent coincide with the roots of 
the equation D £ (z) = 0. 

For 60 > the equation D £ (— A) = 0, for A > 0, has only one solution in (#q/4, 00), then the 
Hamiltonian H £ has only one isolated eigenvalue in (—00, — Oq/A). 

The existence of a pole in the "unphysical" Riemann sheet can be proven by making use of a recursive 
procedure. The equation D £ {z) = can be rearranged as 



iVz - 1 + 1 



2(# + 2iy/z) ' 

To find the solution of the last equation we define the recursive procedure 
Z M = 



+ 2iV^) 4((9 + 2iV^))2 



A: = 0,1,2, 



By using techniques similar to the ones used in the proof of theorem [T] one can see that the sequence 
{z^} is convergent in a ball of radius e 2 around the origin and that the estimate (3.35) holds. 



For 60 = the equation D £ (z) = reads 



e 



0. 



(3.38) 



To find the isolated eigenvalue we set z = —A, A > 0, and rearrange the equation (3.38) as 



4(a/t+a- 1) 



A > 0. 
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Obviously the last equation has only one solution. The recursive procedure 

-4/3 

A (o) 



2 2 / 3 



\(k+l) _ — E ^__ _| fe = 0, 1, 2, 

2v / AW 16A( fc ) 



converges to the solution and can be used to prove the estimate (3.36). 



For 9q < the existence of an isolated isolated eigenvalue and the estimate (3.37) can be proven by 
writing the equation D e (— A) = 0, A > 0, as 



2(Vl + A - 1) 



e 2 



2Vx + \e \ 

and by using the recursive procedure defined by 

a(°) = 



\( k + 1 ) = —Jt + -£ k = 0,1,2,.... 



2\/A~W + |#o| 4(2\/AW + |0 O |)2 



□ 



Let us now analyze the case d = 2. 



Theorem 5. Let d = 2 and assume that c = 0. Then there exists Eq > such that for all < e < Eq: 
the essential spectrum of H e fills the positive real line and is only absolutely continuous, 

a ess (H £ ) = a ac {H e ) = [0, +oo) . (3.39) 

Moreover for e small enough: 

- ifOo > the analytic continuation of the resolvent R e (z) through the real axis from the semi-plane 
Im z > has a simple pole (resonance) in z = E r e and 

4-7TE 2 AlT^ E 2 

Re W = -aWo + ° (e2/l ln£|) ; Im W = ~a%\^ + ° (eVl ln£|2) ' (3 - 40) 

there exists a positive constant C such that the Hamiltonian H e has no isolated eigenvalues in 
(-C,0); 



- if = the Hamiltonian H £ has an isolated eigenvalue in z = E £ and 

E e = -^e 2 |lne| +o(e 2 /\\iie\) ; (3.41) 

- if 0q < the Hamiltonian H e has an isolated eigenvalue in z = E £ and 

E e = -^il + 0{E 2 /\\nE\). (3.42) 

a 1 00 1 

Remark 5. For d = 2, c = and = the Hamiltonian H £ has also one eigenvalue which moves to 
minus infinity as £ goes to zero. 
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Proof. Similarly to the case c < 0, see theorems [T] and [2j for any 6q 6 M. the singularities in z = and 
2 = 1, arising from (—A — z)~ l and (—A + 1 — are absorbed by the terms T Sj h/D £ and Y £ ^2/D £ 
respectively. Then the statement (3.39) is a consequence of the fact that the equation D £ {\) = has 



no solutions for A > and the singularities of the resolvent coincide with the roots of the equation 
D £ (z) = 0. 



For 6q > the equation D £ 



exp [47r( 



l /a)/(9 a 



-A) = 0, for A > 0, has only one solution in (A 0)E3 oo), with A a 
e 2 )] ; this can be seen by setting c 



in equation (3.13). Then the 



Hamiltonian H £ has only one isolated eigenvalue in (— oo, — \ a ,e) 
To prove the existence of the pole E £ in the unphysical Riemann sheet we proceed as we did in theorem 
[Tj We set c = in equation (3.14) and define the sequence 



M = o 

y(*+l) - 



k = 0,1,2,... 



with 



r(z) :-- 



Am 2 



(lnz/(4vr) -l/a-i/4) 



a + (a9 -e 2 )(lnz/(4vr) - l/a-i/4) 



Following what was done in the proof of theorem [l] one can prove that the sequence {z^} converges 



in a ball of radius e 2 and E' £ = z* 00 ) is a solution of the equation D £ (z) = 0. The estimate (3.40) can 
be obtained by computing z^ . 



To find the eigenvalues of H £ for c 
D £ (— A) =0 which reads 



ln(l + A) 

4-7T 



and #o = we set z = —A, A > 0, and study the equation 



e 2 lnA e 2 

—a + — 

Air a 



A > 0. 



This equation can be rearranged as 



ln(l + A) _ 1 
4-7T a 



+ 



e 2 lnA 

47T 



A > 0. 



(3.43) 



The analysis of the r.h.s. and of the l.h.s. of the equation as functions of A is trivial and one can refer 
to what was done in theorem[l]in the study of equation (3.13 ). It comes out that the last equation has 
two solutions. One is in (exp(47r(a/e 2 + l/a)),oo). The eigenvalue of H £ associated to this solution 
goes to minus infinity as e goes to zero (see remark[5]). The other solution is in (0, exp(4-7r/a)) and we 
denote it by A £ . To obtain an estimate of A £ we define the sequence 



a(°) = 

m(l + A( fc+1 ); 
4ir 



1 1 

a + e 2 InAW 
4tt 



k = 0,1,2,.... 



The convergence of the sequence {A^ fc ^} to the fixed point A e can be deduced by examining the 



functions appearing on the two sides of equation (3.43). Such analysis will not be detailed here. 



We limit ourselves to compute the terms which are needed to get the estimate (3.41). By direct 
computation 



A (l) = e 47r/a _ 1 ^/a 



A (2) 



'4vre 2 



exp 



+ 0(e 4 



Aire 2 



A (3) 



^e 2 | lne| + o(e 2 \ lne|) . 



1 /l 



a 2 \a 



ln(e 



Air/ a 



An 



16 



Since the leading order in the expansion does not change in the terms \( k > with k = 4,5,..., the 
estimate (3.41) for the eigenvalue of H £ holds true. 



We conclude the proof of the theorem with the analysis of the case 9q < 0. As usual the eigenvalues 
of H £ are given by the solutions of the equation D £ {—\) = 0, A > 0, which reads 



e 
a 



(a\9 \ + £ 



2 >(1 + A) 

47T 



a+ \9 \ + 



(|e |a+£ 2 )lnA 

47T 



A > 0. 



(3.44) 



Let us denote by /'(A) and / r (A) respectively the l.h.s. and the r.h.s. of the last equation. From the 
analysis of the functions /'(A) and / r (A) one can see that the equation (3.44) has two solutions. The 
function f r has a vertical asymptote in A aj£ = exp [47r(|#o| + a + e 2 /a) /(|#o| a + £ 2 )] and is strictly 
increasing in (0, A aj£ ) and (A ai£ ,oo). Moreover it equals zero for A = 0. The function f l is strictly 
decreasing, it equals £ 2 /a for A = and zero in A = A(°) = — 1 + exp (47T£ 2 /[a(a|#o| + £ 2 )]) • Then one 
solution is in (A aj£ , oo). This solution converges to the eigenvalue of the unperturbed Hamiltonian Hq, 
see the proposition [I] The other solution is in (0, A^). A sequence which converges to this solution is 
defined by A( fc+1 ) = f l {f r (X {k) )), k = 0, 1,2, where we denoted by f the inverse function of / . 
We do not discuss the convergence of the sequence {A^} which can be deduced also by the analysis 
of the plot of the functions f l and f r . The estimate ( |3.42[ ) can be obtained by a direct computation 
of the term A*- 1 ) and noticing that for any k = 1, 2, 3, ... the correction to the term of order e 2 in 
is of order e 2 /\ lne|. □ 

We finally consider the case d = 3. 

Theorem 6. Let d = 3 and assume that c = 0. Then for all e > the essential spectrum of H £ 
fills the positive real line and is only absolutely continuous. Moreover the Hamiltonian H £ has a zero 
energy resonance, i.e., the resolvent R £ {z) has a singularity of order z -1 / 2 in z = 0. 

Remark 6. For d = 3, as opposed to the cases c > and c < 0, when c = the singularity of the 
resolvent Rq{z) in z = does not move from the origin. The perturbation parameter e affects only 
the character of the singularity, turning the embedded eigenvalue into a zero energy resonance. 
Similarly to what happens in d = 1, when #o = 0, there are two solutions of D e (z) = in a neigh- 
borhood of the origin. Both of them are in this case on the negative real axis of the second Riemann 
sheet of y/z. 

We also notice that for 9q < the Hamiltonian H e has an isolated eigenvalue in a neighborhood of 
order e 2 of the point -(47t) 2 /^q. 



Proof. For c = the equation D e (z) = reads 



D £ (Z) 



1 



. #0 r— 



1 



+ 



Wz = 0. 



By a direct analysis one see that the last equation has no solutions on the real positive axes, z = A > 0. 
It is easy to verify that z = is a solution of the last equation for all e > 0. More precisely one can 
see that the function D £ {z) can be expanded around z = as D £ (z) = i(e 2 / 1 {Air) 2 )^fz + 0(|z|). 
Correspondingly, for all e > 0, the resolvent R £ {z) has the following expansion around z = 



R £ (z) = ^ + B £ + 0(\z\ 
v z 



(3.45) 



where A £ is the matrix valued operator with integral kernel 



A £ (x',x) 



( J_JL 

47T \x\ 



1 



i 1 e 



£ \X\ 



i e 
£ I 



-H 1 Ame~\ x \ e-\ x '\ 



x\ \x' 



X 
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According to standard results on the low energy expansion of resolvents of Schrodinger operators in 
dimension three (see, e.g., [IT] and [12J), the presence of a singularity of order 1/2 is the signature 
of a zero energy resonance. In the following we give the explicit form of the resonant state without 



making use of expansion (|3.45|). For this reason, we will not specify any suitable topology in order to 

ning. 

there exists a distributional solution of the equation H £ & r 



give equality (3.45) a rigorous meaning 
for all e ^ 



Let us show that 
Consider the state <3> r given by 



0. 




N 



(- 



\ 



47T 



A 



/ 



where N is an inessential multiplicative constant which we set equal to 1. 
Let us verify that $ r is a zero energy resonance for H £ . The function 



but <t>l i L 2 (M 3 ), 

then <I> r ^ D(H £ ). Nevertheless using formulas given in definition [2j it is possible to compute the 
charges q$ and q\ and the regular parts /q and /{ associated to <3? r . A simple calculation gives 
<?q = — e, q\ = 4tt, /q = and /[ = —1. Since the conditions Qq = ef\ and q\ = —^irfl are satisfied, 
one has that for all * G C$°(R 3 ) Cg°(M 3 ) the equation H £ ^ r ) M - = is satisfied. □ 



3.3 Negative perturbations, c > 0. 

In this section we study the spectral structure of the Hamiltonian H £ in the vicinity of the origin, 



when the parameter c is positive. It is clear form (3.1), (3.2), (3.3), that this choice corresponds to 



perturbations of the Hamiltonian Hq for which the threshold eigenvalue is pushed toward negative 
energies by the perturbative term ce in 9f. 

Theorem 7. Let d = 1,2,3 and assume that c > 0. Then there exists Eq > such that for all 
< e < £ : 

the essential spectrum of H £ fills the positive real line and is only absolutely continuous; 
the Hamiltonian H e has an isolated eigenvalue in z = E £ and 



E e = -ce + 0{e 3/2 ) 



4-7TC 

ce 

27T 



e + 0{e 2 \\ue\) 
+ 0(e 2 ) 



d 
d 
d 



1 

2 
3. 



(3.46) 
(3.47) 
(3.48) 



Remark 7. Notice that when c > the eigenvalue in the upper channel moves toward negative 



energies and is smoothly perturbed by the channel coupling (compare equations (3.46) - (3.48) with 



equations (3.1) - (3.3)). 



For d = 1,2, when #o = 0, the analytic continuation of D e (z) from the semi-plane Im z > through 
the positive real axis has zeros close to the origin. 

As in the case c = 0, see remark |4j in dimension one, one of the roots of D £ (z) = has positive real 
part. The channel coupling produces this resonance as perturbation of the zero energy resonance in 
the lower channel. 

For d = 2 the Hamiltonian H e has also one eigenvalue which moves to minus infinity as e goes to zero. 



The order in e of the remainder in expansions (3.46) and (3.47) are exact for #o = 0. It is possible to 



prove that for 6q ^ the remainders in equations (3.46) and (3.47) are indeed 0(e 



In all dimensions the proof of theorem [7] comes straightly from the analysis of the equation D £ (— A) 
A > 0, and we omit details. 



0, 
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4 Conclusions 



We investigated the spectral properties of model-Hamiltonians describing a quantum particle inter- 
acting with a localized spin via zero range forces. Parameters were adjusted in such a way that the 
unperturbed Hamiltonian showed spectral singularities at the continuum threshold (chosen to be the 
zero energy point). 

All the unperturbed Hamiltonians we considered had a zero energy eigenvalue. In addition, for partic- 
ular values of parameters, a zero energy resonance was also present. Our aim has been to characterize 
the effect of perturbations on the spectral structure around the threshold. 

We defined perturbed Hamiltonians introducing a coupling between the two channels, associated to 
the two possible values of one component of the spin, together with potential-like perturbations of the 
upper channel unperturbed Hamiltonian. 

Direct, sometimes lengthy, calculations bring us to results which agrees with the ones in |10j for 
positive perturbations. In fact the family of Hamiltonians we define have resolvents which are finite 
rank perturbations of the unperturbed resolvent. Roughly, the great part of the considerable work 
done in [10J, in a quite general setting, was to prove asymptotic expansions of the resolvent where 
only finite rank operators appear. This property holds true by construction for all our Hamiltonians 
reducing significantly analytic difficulties and enhancing explicit computability. 

The simplification mentioned above enables us to investigate also purely off-diagonal perturbations, 
where the perturbing term has no explicit bias to move singularities toward larger or lower values of 
energy. Moreover we prove that all the expansions for the singularity coordinates are convergent and 
we give easy recurrent procedure to compute each term in the expansions. 

In particular we want to mention the one and two dimensional pure off-diagonal cases (c = and 
#o = in our notation), when two singularities are present in the unperturbed Hamiltonian spectrum. 
Our results show a peculiar spectral structure of the corresponding Hamiltonian. As expected, no 
continuity of the spectral properties in parameter space is observed in this particular point. 
In this paper we analyzed only multi-channel Hamiltonians. As we mentioned in the introduction, one- 
channel Schrodinger operators describing a quantum particle interacting with many point scattering 
centers show very rich spectral structures and can suitably approximate Hamiltonians with any kind of 
smooth potentials. Moreover their resolvents are finite rank perturbations of the Laplacian resolvent 
for any (finite) number of scattering centers. In our opinion such kind of Hamiltonians are good 
candidates to examine spectral properties of a vast class of Schrodinger operators. 
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